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ABSTRACT 

We examine the evolution of residual gas within tidally-limited dwarf galaxies and globular 
clusters. In systems where the gas sound speed exceeds about 10% of the central velocity dis- 
persion, hydrostatic equilibrium causes the gas to have significant density at the tidal radius. In 
contrast to the stellar component, gas particles relax rapidly in response to mass loss. Gas lost 
through the Lagrangian points is quickly replenished, leading to continuous outflow, analogous to 
the Jeans' escape of planetary atmospheres. In such systems, the gas may be lost on timescales 
as short as a few times the sound crossing time of the system. In colder systems, by contrast, the 
density contrast between the core and the tidal radius is much larger, greatly reducing the mass 
loss rate, and allowing the system to retain its gas for over a Hubble time. Using this criterion as 
a guide, we present the results of numerical simulations to show that either in the proximity of 
large galaxies or in clusters of galaxies, residual gas in typical dwarf spheroidal galaxies (dSphs) 
and globular clusters would be efficiently lost after it is ionized and heated, either by external 
or internal sources. The tidally removed gas shall follow an orbit close to that of the original 
host system, forming an extended stream of ionized, gaseous debris. Such tidal mass loss severely 
limits the ability of dSphs and globular clusters to continuously form stars. The ordinary gas 
content in many dwarf galaxies is fully ionized during high red-shift epochs. In extreme cases, 
therefore, star formation may be completely prevented by early segregation of dark and ordinary 
matter, leading to the formation of starless, dark-matter concentrations as the remnants of dwarf 
galaxy mergers in the CDM model of galaxy formation. In the field, where the dwarf galaxies are 
tidally isolated, and in the center of clusters of galaxies, where the tidal potential is compressive, 
ionized gas may be retained by dwarf galaxies, even though its sound speed may be comparable 
to or even exceed the velocity dispersion. In such situations, dwarf galaxies may retain their gas 
content until they merge with much larger galaxies. These processes may help to explain some 
observed differences among dwarf galaxy types, as well as observations of the haloes of massive 
galaxies. 

Subject headings: galaxies: dwarf — galaxies: evolution — clusters: globular — clusters: open 
— hydrodynamics — methods: numerical 



1. Introduction 

In the widely adopted cold dark matter (CDM) cosmological models, small objects are the first to form, 
and larger galaxies form via the merger of dwarf systems (e.g. White & Rees 1978; Blumenthal et al. 1984; 
Cole et al. 1994; Klypin, Nolthenius, & Primack 1997; Navarro, Frenk & White 1997). The evolution 
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of the smaller building blocks is therefore of great interest for the evolution of larger systems. Because 
they formed prior to the emergence of normal galaxies, these subsystems are theoretically expected to have 
greater density than the larger systems, and therefore they are likely to retain their integrity in the halo 
of large host galaxies, which would imply the existence of many more surviving systems than are observed 
today (Moore et al. 1999). The lack of small stellar systems around larger host galaxies has been referred 
to as the "missing satellite" problem (Klypin et al. 1999). 

One possible resolution of the problem is that gaseous baryonic matter may have been ejected from the 
dwarf systems due to massive star formation, which in turn leads to photoionization-driven outflows, stellar 
winds, and supernovae (Tenorio-Tagle et al. 1986; Efstathiou 1992; Lin & Murray 1994; Quinn, Katz, & 
Efstathiou 1996; Kepner, Babul, & Spergel 1997; Navarro & Stcinmctz 1997; Weinberg, Hcrnquist & Katz 
1997; Barkana & Loeb 1999; Benson et al. 2002; Dong, Murray, & Lin 2003; Dekel & Silk 1986; Mac 
Low & McCray 1988; Madau, Ferrara, & Rees 2001; Mori, Ferrara, & Madau 2002). Once segregated, the 
ordinary matter continues to its dissipational infall, while the dark matter components form a lumpy halo 
in the potential of the large galaxy (Weil, Eke & Efstathiou 1998; Sommcr-Larscn, Gelato, & Vcdcl 1999). 
If the mass loss is driven by bursts of star formation, a small amount of residual ordinary matter, which has 
already been converted into stars, would remain with the clumpy dark matter subsystems. 

In the halo of the Galaxy, several dwarf spheroidal galaxies (dSphs) are found in the form of diffuse 
patches of stellar light (Hodge 1971) embedded in dark matter haloes (Aaronson 1983; Faber & Lin 1983). 

It is natural to identify these with the remnant subsystems, i.e. the theoretical building blocks of large host 
galaxies in the CDM scenario (Mateo 1996; Klypin et al. 1999; Moore et al. 1999), because they have 
similar total mass (> 1O^~^M0). There are, however, significant differences. Both the small (few km s~^) 
velocity dispersions, and the inferred dynamical density of the observed dwarf galaxies arc smaller than those 
postulated in the CDM models. Also, the spread in the observed amount of luminous matter (and hence the 
mass to light ratio) in the dSphs is much larger than that in their total mass (Mateo 1998). 

Some of the above differences may be due to evolution. The dSphs are found to have a wide range 
of star-formation and enrichment histories, with few clear trends (Mateo 1998; Grebel 2001; Harbeck et 
al. 2001). All contain populations of old stars of varying sizes, and ages exceeding 10 Gyr. In addition, 
some show histories of continuous star formation, many others show evidence of bursts of star formation 
separated widely in time, and a few show little or no evidence of subsequent star formation. The observed 
stellar population and the modest mass-to-light ratio of the dSphs are consistent with the hypothesis that 
a significant fraction of the original baryonic matter may have been ejected from them and other dwarf 
galaxies as a consequence of earlier epochs of star formation (Gerola, Carnevali, & Salpeter 1983; van Zee, 
Skillman, & Salzer 1998). 

In addition to dSphs, the families of low-mass, gas-poor dwarf stellar systems also include dwarf elliptical 
galaxies, and globular clusters. These groups of objects are quite distinct from each other (Kormendy 1985). 
Some of their differences may, however, result from post-formation evolution processes, such as mass loss 
either prior to or following star formation (sec discussion above), or segregation of baryonic and dark matter 
during star formation (Murray & Lin 1989). The globular clusters are the low-mass extreme of these families 
of dwarf stellar systems. Their velocity dispersions are generally more than twice that of the dSphs within 
the Local Group of galaxies, but less than that of typical dwarf elliptical galaxies within nearby clusters 
of galaxies. In contrast to the dwarf spheroidal galaxies, all but one or two globular clusters of the Milky 
Way show evidence for only a single, old stellar population. As in those dwarf galaxies which show little 
evidence for more than a single generation of star formation, any gas within globular clusters that remained 
after the initial generation of stars formed must have been lost from within them. Nevertheless, the multiple 
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populations of stars in at least one globular cluster [u Cen) and several dwarf galaxies indicate that the star 
formation process may not always be an effective mechanism for decoupling the ordinary and dark matter 
in these loosely bound stellar systems. 

Collectively, the population of the observed dwarf stellar systems falls far below that predicted in the 
CDM model, and so they alone cannot be the resolution of the "missing satellite" problem. Gravitational 
lensing observations, however, find preliminary indications of starless dark matter clumps in the halos of 
nearby galaxies (Mctcalf 2002; Mctcalf & Zhao 2002). If such dark entities dominate the halo structure, 
then an efficient and severe segregation of ordinary and dark matter, through processes in addition to the 
violent impact of star formation, is needed. A more quiescent process may therefore also play a significant 
role in determining the ability of dwarf stellar systems to retain their gas to form multiple generations of 
stars (Lin & Murray 1998). We shall address the cause of multiple generations of stars in dwarf spheroidal 
galaxies in a separate paper. 

Several gas dynamical processes may lead to gas removal from stellar systems with shallow potentials. 
Prior to star formation, gas may be removed by either ram-pressure or Kelvin- Helmholtz instability resulting 
from motion through the halo gas of a larger system (Lin & Fabcr 1983; Murray et al. 1993). The efficiency 
of gas removal due to either process depends upon the depth of the gravitational potential, and the density of 
ambient gas. While it is possible that the dwarf galaxies and globular clusters seen today were able to retain 
sufficient gas to form stars within them because they have deeper potentials than the dark matter-dominated 
"missing satellites," the amount of ambient gas in the halo of the large, host galaxies is highly uncertain. 

In the halo of the Milky Way and other galaxies, both satellite dwarf galaxies and star clusters are 
subject to tidal forces. For dSphs, this process occurs following the formation of a nearby, massive galaxy, 
whereas star clusters form within the potential of a massive system. The tidal forces increase rapidly in 
significance for smaller orbits around the host galaxy, and can strongly perturb any ionized gas within the 
satellite stellar systems. The resulting situation is somewhat analogous to that of interacting binary stars 
(Pringle 1985), in that the ionized gas in the satellite systems always relaxes, and expands to fill their Roche 
Lobes. A fraction of the gas continually streams away from the satellite systems through their inner and 
outer Lagrange points. If the rate of gas loss is substantial, it may limit the ability of dwarf systems to 
retain gas for subsequent generations of star formation. 

In this work, we examine the significance of mass loss from tidally- limited dwarf systems. We begin 
in Section 2 with an analytic estimate of the rate of mass loss. In Section 3, we examine the results of 
simple one-dimensional models, and continue with three-dimensional models in Section 4. Our results are 
summarized and consequences for the evolution of dwarf systems are drawn in Section 5. 



The distance of the LI and L2 points from the center of the dwarf satellite are given approximately by 
the Roche radius. 



2. Analytic Determination of Mass Loss Rates 




(2-1) 



for a point-mass parent-galaxy potential (King 1962), and 




(2-2) 
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for a logarithmic parent-galaxy potential (Oh & Lin 1992). In the above expression, is the total mass 
(gas, stars, and dark matter) of the satellite system. The parent galaxy has mass Mq inside the orbital 
semimajor axis a of the satellite system. Although the concept of the tidal radius for a satellite with an 
orbital eccentricity e is uncertain and controversial (Innanen 1979; Allen & Richstone 1988), it provides a 
useful fiducial length scale (Oh & Lin 1992). The above results may be parameterized as 



pMd + Ma 

where /? is a function of e. For circular orbits around a point mass potential, /? = 3, whereas around an 
isothermal potential, (3 = 2. 

The rate of mass loss through the Lagrange points is related to the Roche radius and the sound speed 
of the gas, c^. A full derivation in the context of a low-mass planet orbiting around a main sequence star is 
given in Gu, Lin & Bodenheimer (2002). For the total loss rate of gas through the LI and L2 points, the 
results can be approximated as 

Mg = AwpiRlcJ, (2-4) 

where pi is the gas density at the Roche radius. The geometrical correction factor, /, can be approximated 
by expanding the potential near the Lagrangian point (see 4-3). 

If the gas within the dwarf system is isothermal, and the potential of the system is approximated with 
a Plummer potential (Binney & Tremaine 1987), then the gas density varies as 



Pg{r) =Poexp 

while 



- i^r) - <^o) 



-GMa 



(2-5) 



{R^ + r2) 

where (/)o = 4>{r = 0) = —GM^/ Rc- In the above, Rc is the core radius, defining the radius over which the 
density of gravitating matter is approximately constant. In this potential, the ratio of the gas density at Rr 
to that at the center (po) is 



Po 

where xr = Rr/Rc, and 



exp 



(1+4)1/2 



(2-7) 



Even for relatively concentrated systems where xr 1, the gas density retains a finite value at large radius, 
albeit p{Rr)/po can be very small for cold gas, with large values of a. For large values of p{Rr), self-gravity 
of the gaseous envelope may become important, and the gas shall become a self-gravitating isothermal sphere, 
resulting in a more rapid decline in density than assumed here. 

In response to mass loss, the gas adjusts on a sound crossing time in an attempt to re-establish a 

density distribution appropriate to hydrostatic equilibrium. Any loss from the system at large radii is 
therefore continually replenished. This is in contrast to the stars, which diffuse on the much-longer two- 
body relax;ation time scale. Thus, the mass loss rate of the gas, Mg remains approximately constant prior 
to any significant reduction in Mg. 
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The total mass of isothermal gas within the dwarf system is given by 



Jo 

where x = r/Rc. For a Plummer potential, 

Mg = AnpoRl f 
Jo 



exp 
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x'^dx. 



(2-9) 



(2-10) 



Note that, for sufficiently large x^, the dwarf galaxy's potential can contain an arbitrarily large mass, even 
in the hot limit where a <1. But for modest values of either Xn or the mass of gas within the potential 
is limited. 

Using Equations (2-2), (2-4), and (2-10), and assuming the dwarf satellite to be on a circular orbit of 
period 

, .-3 

P = 2n\ 



the mass loss time scale becomes 

where 

r {xr, a) = exp < a 1 
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(2-11) 
(2-12) 
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gives the number of orbits of the dwarf satellite around the parent system before significant mass loss occurs. 
In the limit of cold gas, a :s> 1 and 
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so that over a Hubble time scale, no significant mass loss occurs. But, for a < ai, where 

1 



ai 



■ exp 



ai 



(1 + 4)'^', 



(27r/3)V24/V, 



(2-14) 



(2-15) 



a significant amount of gas is lost within an orbital period. 



In Figure 1, we show the variation of F {xR,a) with a (panel A) and with a' (panel B) for three values 
of xr, selected to match the values used in the one- and three-dimensional dwarf galaxy models discussed 
in the next two sections. For these calculations, we choose /3 = 2, corresponding to a circular orbit in an 
isothermal host-galaxy potential. The quantity a' is defined as 



a' = a 



(2-16) 



such that a' corresponds to the Mach number of the escape speed from the center to the Roche radius. 

As can be seen from panel A of Figure 1, and as would be expected from Equations 2-12 and 2-13, the 
time scale for loss of gas may become extremely short, less than the orbital period, if a < ai ~ 10 — 15. The 
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curves of T{xji, a) vs. a for xji = 1, 3, and 9 cross each other for large a. We note, however, from panel B 
of Figure 1, that for a > 1, T{xn, a) is a more sensitive function of a' than of a. The critical condition for 
rapid mass loss, r{xR, a) = 1, occurs when a' = a[ o±5 — 10, depending on the value of xr. 

The observed velocity dispersions of dwarf stellar systems, which are indicative of the potential depths, 
are similar to the sound speed of ionized gas (Mateo 1998; Pryor & Meylan 1993). Once the gas in these 
systems is ionized, either by internal or external sources, it shall have a' ~ 1. In this limit, the results in 
Figure IB imply that the ionized gas may be rapidly lost by a continuous, quiescent outflow, thus preventing 
the formation of subsequent generations of stars. Subsequent stellar generations would form only if gas could 
be accreted, a possibility which shall be examined in a subsequent paper. 



3. One-Dimensional Models 

In the analytic model presented above, the gas is assumed to bo in quasi-hydrostatic equilibrium at all 
times. Yet, mass loss leads to outflow and a reduction in the residual gas mass which should be treated 
with a time-dependent dynamical calculation. We examine the time evolution of the gas in dwarf galaxies 
and star clusters using numerical simulations. In order to compare with the analytic results above, we first 
carry out a series of one-dimensional dynamical simulations. Using these models, we show the importance 
of ionization and gas temperature on the rate of tidally induced outflow. 

The one-dimensional scheme used was first developed to study star formation and feedback effects in 
dwarf galaxies (Dong, Murray, & Lin 2003). It is a ID Lagrangian code, into which radiative transfer and 
photoionization are incorporated using a Stromgrcn shell model. It thus allows us to explore the effects of 
heating and photoionization by both internal and external UV sources. With a Lagrangian scheme, we also 
do not need to specify artificial conditions near the boundaries of the computational domain. Gas is allowed 
to freely flow out of the region. 

For the ID models, the gas experiences the gravity of the dark matter of the dwarf galaxy, as well as 
an external tidal potential. The dark matter potential of the dwarf galaxy is taken to follow the form found 
by Burkcrt (1995), which has been tested to be a slightly better fit to dwarf spheroidal galaxies than the 
simpler Plummer potential, used in the previous section. The corresponding dark matter density distribution 
is given by 

3 

PDM = 7 — ; TT^— — 2T- V"5-lj 

(r + ro)(r^ + r^) 

The central density and velocity dispersion of the dwarf galaxies are scaled with a structural length scale ro 
such that 

Pe = 4.5 X 10-2(ro/kpc)-2/=^MQpc-^ (3-2) 
a = 17.7(ro/kpc)2/3km s'^ (3-3) 

The gravitational potential due to the dark matter within the dwarf galaxy is determined from its density 
using Poisson's equation. From equations (2-5) and (2-10), we can determine the density distribution of an 

isothermal gas in this potential, including its central value po. In order to determine the extent of the dwarf 
galaxy, we also impose an external tidal potential, chosen to follow the form 



<l>tide = <t>0 
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whcre (jjQ < 0. The exact value of A is found be to unimportant; in the following calculations, we choose 
A = 1. We specify the tidal radius Rr to be the location where the gradient of the total potential vanishes. 
One can show that Rr/to is related to pc and 0o as 




Rr 
ro 



— arctan 



Rr^ 



(3-5) 



The initial distribution of gas is set up in hydrodynamic equilibrium under the dark matter potential 
alone. The total mass of gas, Mg = qMjjM, where the ratio of baryonic to dark matter q is assumed to be 
0.1 in most of our calculations (Mateo 1998). At the onset of the calculations, both the external potential 
and extragalactic UV radiation are switched on instantly . The residual gas quickly settles into a quasi- 
equilibrium after a brief and minor adjustment in the outermost region. Because the gas density is tenuous 
in this region, this adjustment does not affect the subsequent evolution of the gas. 

We assume the gas to be metal-frec, so that in the absence of photoionization it would cool to a 
temperature T ~ 10'^ K (Kang et al. 1990; Lin & Murray 1992). The temperature of the photoionized 
gas is taken to be 15,000 K, appropriate for low-metallicity, photoionized gas. The external UV flux Ji, w 
10~^^ ergs s~^ cm~^ ster~^ Hz~^ at the Lyman limit at low redshift, and 50 times more intense at redshift 
z = 3 (cf. Haardt & Madau 1996), is sufficient to completely ionize most dSphs, as confirmed by our models. 
More massive dwarf elliptical galaxies (dEs), are, however, partially self-shielded from the external UV, 
and so they contain a region of neutral residual gas surrounded by ionized gas. Because we are, in this 
section, primarily interested in the efficiency of photoionization for dEs, the model parameters used here 
are appropriate to those systems. In the next section, we present the results of some 3-D models for the 
evolution of fully ionized gas within dSphs. 

The properties of the one-dimensional models are shown in Table 1. For each model, we list M^, ro, the 
escape speed from the central potential {a, equivalent to GMd/Rc in the previous section), a' for the ionized 
gas, J^, the initial baryonic mass fraction (Mgi/Md), the ratio of the final mass of neutral gas to the initial 
total gas mass {Mh / /Mgi), the ratio of the initial radius of the sphere of neutral gas to the Roche radius 
{Rii/Rn), P, the computed mass loss timescale for gas {tioss), the ratio of tioss to P, and the computed 
value of r. 

In the standard model (lD-1), ro = 1 kpc, such that cr ~ 18 km s~^. In this and all other models 
presented here, wo set Rr = 3ro, so that the dwarf galaxy's total mass Md 5Mo where Mq is the total 
mass within ro. In Model lD-1, the galaxy's mass (including both dark and gas components) within ro and 
Rr are, respectively, Mo = 7.2 x IO'^Mq and Md 3.66 x 10^ Mq. 

If the dwarf galaxy has a circular orbit in an isothermal potential of the background host galaxy, with 
circular speed Vc = 210 km s~^, and mass 

Mg = 5 X lO"(a/5Okpc)M0, (3-6) 

the orbital radius a is given by Equation (2-3). The orbital period, P = 1.49 (a/50 kpc) Gy and (3 = 2. 
For Models lD-1, lD-2, lD-3, and lD-4, we adopt ro = 1 kpc so that a = 129.1 kpc and P = 3.85 Gy. For 
Model lD-5, we choose ro = 1.35 kpc, such that a = 182.6 kpc and P = 5.14 Gyr. 

The evolution of the mass fractions of gas within the tidal radius for the one-dimensional models are 
shown in Figures 2 and 3. As can be seen in Figure 2, half of the gas mass is lost from Model lD-1 within 
350 Myr, in spite of the fact that Cg — 0.6a. This result is consistent with our discussion in the previous 
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section where we showed that a significant fraction of the ionized gas is lost within one galactic orbital period 
for a = (cr/cg)^ < ai ~ 10. Such rapid outflow of gas would prevent extended and protracted star formation. 
The ionized gas at the outskirts of the galaxy flows beyond Rr with the sound speed (~ 10 km s~^) of the 
ionized gas. 

In Model lD-1, = 4 x 10~^^ ergs s~^ cm~^ ster~^ Hz~^, appropriate for redshift z = 1. In this 
model, a small amount, M^f, of the residual gas remains in a neutral state at the center of the galaxy, 
surrounded by ionized gas at an interface of radius i?„. In Figure 4 are shown snapshots of the densities 
of the neutral and ionized gas, p„ and pi for the first four one-dimensional models, as a function of radius. 
These results show that, analogously to the Stromgren sphere in an HII region, the ionization front contracts 
slowly from the outer tenuous to the inner, denser regions of the galax;y, while the ionized gas lost from the 
outer halo is being continually replenished. 

When a quasi-static equilibrium is established, in which the external UV flux is balanced by rccombi- 

1/2 

nation in the outer ionized envelope, the magnitude of pi oc J^' (Figure 4), while Cg is independent of J^. 

1 /2 

We thus expect M oc Ju ■ This conjecture is confirmed by Models lD-2 and lD-3 which have identical 
structural parameters as Model lD-1, but reduced J,y. In these cases, both the mass fraction and radius of 
the neutral sphere are larger than the standard model. Similar to the standard model, a ~ 1 < ai in the 
ionized envelope and the ionized gas flows from the ionization front at i?„ to Rr in a fraction of the sound 
crossing time. The density of ionized gas at Rr is, however, reduced relative to the standard model, due 
to the larger masses of cold, neutral gas, with its steeper radial density falloff. The result is a reduction 
in Md relative to the standard model. The actual dependence of J^, upon redshift remains uncertain, and 
the values listed in Table 1 are comparable to the values estimated for the present epoch by some authors 
(Haardt & Madau 1996; Bechtold et al. 1987; Ostriker & Ikeuchi 1983). 

In Model lD-4, we increase the ratio, q, of the baryonic-to-dark matter mass by a factor of two from 

the standard model. The increased density of baryonic matter enhances self-shielding, and the gas becomes 
mostly neutral. The gas cools and settles to the center of the galaxy, with a smaller The ionization rate 
across the front is balanced by rapid recombination, making photo-evaporation less efficient. Consequently, 
the value of pi decreases, and the mass loss rate is reduced (Figure 3). In Model lD-5, we increase tq by a 
factor of 1.35 but maintain the baryonic fraction at 0.1. In accordance with the Burkert (1995) potential, 
Md in Model lD-5 is twice that in Model lD-1, though the initial gas mass Mgi is the same as Model lD-4. 
With these parameters, gas in the center of the dwarf galaxy is less dense, but the total baryonic mass is 
greater than in the standard model. Once again, the larger gas reservoir enables a greater fraction of Mgi to 
become self- shielded, neutral, cool, and settle to the center of the galaxy, again reducing the mass loss rate 
relative to the standard model. 

We have also computed a case which has identical physical parameters as Model lD-1, but the galaxy 
is set in isolation (infinite tidal radius). In this case, the gas adjusts to both ionization and hydrostatic 
equilibrium with no gas attaining finite outflow velocity. As indicated above, even for the early epochs 
where z and are large and the gas is fully ionized and hot, these static equilibria can be attained from 
Equations (2-5) and (2-10) with a sufficiently large xr or equivalently Rr. The implications of these results 
will be discussed in § 5. 
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4. Three-Dimensional Models 

The one-dimensional models provide the basic description for the underlying ionization, relaxation, and 

hydrodynamics of the residual gas which regulate the efficiency of its tidally induced outflow. But in reality, 
the tidal forcing of the host galaxy is non-spherically symmetric about the center of the dwarf galaxy. In order 
to consider a more realistic tidal potential, we examine the results of some three-dimensional simulations in 
this section. 



4.1. Numerical Method 

The numerical code wc use is Cosmos. It is a massively parallel, multidimensional, radiation-chcmo- 
hydrodynamics code developed at Lawrence Livermore National Laboratory. The code, and tests of the code 
are described in Anninos & Fragile (2003), and Anninos, Fragile & Murray (2003). It has also been used to 
study the evolution of supernova-enriched material in dwarf galaxies by Fragile et al. (2003) . The reader is 
referred to those papers for more details on Cosmos. We discuss here the settings used for the current work. 

Wc wish to examine the evolution of cithcT cool or heated gas within dwarf galaxies and massive star 
clusters. The one-dimensional models discussed above include the effects of heating, cooling, and ionization 
of the gas and show that the gas dynamics is primarily determined by its temperature and more indirectly 
by the state of ionization. For the three-dimensional models, the inclusion of ionization and recombination 
of the gas would only reduce the controlled nature of the models, without adding new information. Thus, we 
choose to set the gas to be initially isothermal, at a determined temperature. In order to approximate the 
effects of heating and cooling, which has the tendency to maintain the gas at a nearly constant temperature, 
we assume that the gas evolves adiabatically, with exponent 7 = 1.1. We choose two sets of models. For 
the hot models, we set the gas temperature to be 8,000 K and for the cold model we assume the gas is in 
an atomic state with a temperature of 1,000 K. These isothermal states correspond to either early or recent 
epochs, when is either adequate or insufficient to significantly ionize the residual gas. The results of the 
1-D models in the previous section indicate that, for dSphs, the entire galaxy is likely to be ionized by a 
modest external UV flux. Thus an isothermal temperature is a reasonable approximation for the hot models. 
When the gas attains a sufficiently large density, it becomes self-shielded inside these galaxies and clusters. 
The isothermal approximation is also sufficient for this cool phase. 

Self-gravity of the gas is not included in the models. This approximation should be adequate for dwarf 

galaxies, which arc dark-matter rich, and so their gravity is dominated by that of the dark matter. For a 
globular cluster, which has undergone efficient star formation, the potential is dominated by the stars, and 
so, again, the potential of the gas is not important to the evolution. The models therefore include a static 
background potential, within which the gas evolves. In order to compare the numerical results with our 
analytic calculations, the background is modeled as a Plummer potential, as expressed in Equation (2-6), 
which has a form similar to, but simpler than, those found to be followed by dwarf galaxies (Burkert 1995), 
and which has also been used to model globular clusters (Plummer 1911). 

To the potential of the dwarf stellar system is added that of a massive parent galaxy. The total density 
of the massive parent galaxy is assumed to follow the same isothermal form as above, with a total mass of 
5 X 10^^ M0 and a circular speed Vc = 210 km s~^ at some typical distance 50 kpc from the Galactic center. 
Smaller Galactic radii arc used for the globular cluster models. 

Throughout the numerical simulations, our models remain centered upon the dwarf systems, ie. we are 
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in a non-inertial frame of reference. Coriolis terms must therefore be added to the potential, which then has 
the form ^ 

(j>{x, y, z) = (j>d{x, y, z) + (pcix, y, z) - -fJ^ [x^ + y^) , (4-1) 

where <f)d,G are the conservative gravitational potentials of, respectively, the dwarf stellar system and parent 
galaxy, is the angular velocity of the rotating frame (the system is assumed to be in a circular orbit), the 
2:-axis is the axis of rotation, and the a;-axis points towards the center of the parent galaxy. This form of the 
effective potential is analogous to that used for close binary stars (Pringle 1985). 

The form of the potential is shown in Figure 5, which shows a surface plot of (/)(a;, y, 0) for Model 3D-ld. 
Positions on the x and y axes, and the potential are all in code units. The LI and L2 points are located 
along the a;-axis, at a; = ±0;^, where Xt = 0.89 in code units. Due to the saddle-shape of the potential at LI 
and L2, mass loss may occur only near the Lagrange points. At an orbital radius a where <C V^a/G 



Xt = ± 



GMd ^ 



1/2 

Rc (4-2) 



which reduces to xt ~ (GMd/2n'^)^^^ in the limit that Rc <C (GAfd/20^)^/^. In comparison with our 
analytic, spherically symmetric approximation, xt is the analogous quantity to Rr, and is only slightly 
modified from the expression for the Roche radius of a point-mass potential for the dwarf galaxy. 

In the vicinity of Xt, the gas density is approximately constant because the gradient of ^ vanishes, such 
that the density scale height of the gas diverges at Xt- The escaping gas streams out of the dwarf galaxy 
through the saddle surfaces, which are bounded by the local potential maxima in which d(j)/dx = 0. This 
surface lies approximately on the y — z plane, centered on a; = Xt- For Xt > GMd/G^, we can estimate 
the magnitude of / (cf. Equation 2-4) under the assumption that gas diffuses to the outflow channels in a 
hydrostatic manner, such that 



1 r°° r 

47ra;? J_, 



f ~ 7Z72 / / exp 



4>{xt,y,z) - 0(a;t, 0,0)' 



•^•"^r-(-¥)-r(4)-(i)"'(si)^?- 



dydz 



(4-3) 



Although hydrostatic equilibrium is approximately maintained within the tidal radius, the outflow via the 
channels prevents the gas from establishing a sustained equilibrium, and that the magnitude of / may be 
quite different from the above estimate. Nevertheless, the above approximation suggests that / is unlikely to 
be orders of magnitudes smaller than unity for modest to large values of a' . For a < 1, / in Equation (2-4) 
cannot exceed unity. 

In principle, non-inertial terms should also be added to the equations of motion of the gas. These would, 
however, strongly affect only that gas which was either at or beyond the tidal radius. Because our goal in this 
work is to examine the rate of mass loss, and not the subsequent evolution of gas lost to the dwarf system, 
omission of the non-inertial terms shall not affect our conclusions signiflcantly. Our models are generally 
integrated for a small fraction of an orbit, and so their affect would be small in any event. 

The gas density is initially set up such that 



p{x, y, z) = po exp [-(/> {x, y, z) /c^] , 



(4-4) 
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where is the isothermal sound speed of the gas. For a spherically-symmetric potential, such an initial 
distribution of isothermal gas would describe hydrostatic equilibrium. For the non-spherical potential de- 
scribed above, however, the density of the gas near R = xt evolves significantly at early times. In order 
to prevent this from causing us to overestimate the mass loss rates, we calculate the mass of gas contained 

within R < xt/2, which wc denote as M2(<). That quantity shows little affect from the behavior at larger 
radii. In addition, we avoid influences from early behavior by fitting the value of M2{t) at late times. 

Cosmos is written for Cartesian coordinates, and so we cannot take advantage of the symmetry of the 
potential to run the problem in two dimensions. In order to reduce computation time, the problems are run 
as quadrants. Typical resolutions are 200x100x100. We have examined the eff'ects of changing resolutions 
in some of the models, and find that doubling the resolution changed the time scales for mass loss by only 
10%. 



4.2. Model Results 

The models which we have run in order to examine the evolution of gas within dwarf galaxies are shown 
in Table 2. On the table are listed, for each model, the total mass of dark matter in the dwarf. Mdm, the 
core radius of the dark matter potential, i?c, the depth of the potential, (pQ^^ (measured from the center to 
the galaxy to xt along the y-axis, i.e. (j)o [l — (1 -|- a^/j)""^^^] , the orbital radius of the dwarf galaxy around 
the parent system, a, the value of Xt (i.e. Rr), the initial gas temperature, Tq, a ', xr{= Xt/Rc), P, 
the measured mass loss time scale, tioss, the ratio PT(xb, a')/Tioss, and the equivalent / as derived from 
Equation (4-3). In the derivation of F, we set /3 = 2, as in Figure 1. The galaxy parameters are chosen to 
be representative of small dwarf galaxies (Mateo 1998), and the variations of the parameters examine the 
effects of mass, orbital radius, and gas temperature upon the mass loss rate. The high temperatures here 
and below arc chosen to be representative of either warm neutral, or warm ionized gas with metal cooling. 
Throughout the simulation, the gas temperature remains nearly constant. 

The values of M2(t)/M2(0) are shown in Figure 6. In deriving the Vcilues of tiqss 

quoted in Table 2, the 

curves have been fitted to the form 

M{t) = M(G) exp i-t/Tioss) . (4-5) 

As can be seen from Figure 6 and below, the curves are not always well represented by exponentials. It 
is apparent, however, that the derived values of tioss represent well the rates at which gas is lost from the 
systems. The time scales for mass loss are comparable to that seen in the one-dimensional case above. In 
comparison with the analytic results in §2, we find from Equation (2-7) that / ~ 0.2 — 1.2 which agrees poorly 
with that inferred from Equation (2-12). As indicated above, the departure from hydrostatic equilibrium 
may be the cause of this disagreement. In Figure 7, we plot the mass fiow field (pv) for Model 3D-4d at 
680 Myr. The figure shows a slice in x and t/, at z = 0. The greatest values of jpvj occur in the complex flow 
at the center of the galaxy. Away from the center, the vectors of pv are almost parallel to the a;-axis. The 
semi-circle in Figure 7 has radius Xt, and indicates the size of the Roche lobe in the spherical approximation 
(the actual potential used in the model is highly aspherical, as indicated in Figure 5). Clearly, only a fraction 
of the Roche surface is open to outflow. The complex flow within the core is timc^-dcpendcnt, but the flow 
outside of the core is not observed to vary significantly over time. The largest velocities present, near the 
tidal radius, are 10 km s~^. 

Also apparent from the differences of tioss among the models is the extreme sensitivity to a (cf. Equa- 
tions 2-12 and 2-13). The critical value of a ' ^ 6, half the value derived in the analj^ic approximation. For 
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a < ai, a significant fraction of the original gas in the dwarf galaxy is lost within an orbital period. The 
ionized gas in both the lower-mass galaxy and the higher-mass system at lower orbital radius (Models 2 and 
4) is largely lost within 1 Gyr, which would preclude extended epochs of star formation. 

The globular cluster models are listed in Table 3, which lists the same quantities as in Table 2. The 
masses and core radii are representative of globular clusters (Pryor & Meylan 1993). The orbital radii are 

either just outside of the Solar circle, or well within it (most globular cluster orbit closer to the center of the 
Galaxy than the Sun). The evolution of M2(t)/M2(0) for the models is shown in Figure 8. 

As was the case for the dwarf galaxies, ionized gas is lost quickly from globular clusters, on time scales 
< 100 Myr. The critical value a' ~ 6 as for the dwarf galaxy. Again, the rapid mass loss induced by 
tides precludes extended epochs of star formation in these systems. Note that for a' ~ 1, the theoretically 
determined / is much larger than that estimated from (2-12). In this limit, the assumption of quasi- 
equilibrium breaks down near the Lagrangian point. Nevertheless, the dependence of tioss on the magnitude 
of a' is well established. 



5. Discussion and Consequences for the Evolution of Dwarf Systems 

The essence of our results is contained in Tables 1, 2, and 3. The results of the calculations in §3 indicate 
that gas in typical dwarf spheroidal galaxies (dSphs) can be fully ionized, even by a modest background UV 
flux. But for the more extended and dwarf elliptical galaxies, a fraction of the residual gas can remain neutral 
at their central regions, especially during more recent epochs. Even for these galaxies, mass loss occurs if 
the dwarf galaxies are captured as satellites of larger galaxies or clusters of galaxies. 

As can be seen in the tables, and expected from § 2, the time scale for mass loss is extremely sensitive 
to the ratio a = (/)o/c^. The magnitude of (/)o can be directly measured from the central velocity dispersion. 
If a » ai ~ 6 — 10, then the gas is tightly held within the potential of the dwarf stellar system. The density 
contrast between the center of the galaxy and the tidal radius, given by Equation (2-7), p{xt, 0, 0)/p{0, 0, 0) ^ 
exp(— a) ^ 1. In this limit, the time scale for mass loss, tioss ^ ^d, where is the dynamical time scale of 
the dwarf system. The reverse is true when a = (j)o/cl < ai. 

For many dwarf galaxies, the velocity dispersion is comparable to the sound speed of ionized gas, and 
so, if the gas is ionized, the system should have a « 1. Once a < ai, our results indicate that the gas 
may be lost to the dwarf system on a time scale no longer than several times r^, which is shorter than 
both the orbital period, P, and the time scale on which subsequent generations of stars would be expected 
to form (Murray & Lin 1993; Lin & Murray 1994). Multiple possible sources of ionizing radiation exist, 
including extragalactic UV, UV from a nearby massive galaxy, or UV from stars formed within the dwarf 
galaxy (Weinberg, Hernquist & Katz 1997; Kepner, Babul, & Spergel 1997; Dong, Murray, & Lin 2003). 
The evolution of dwarf galaxies therefore depends upon their proximity to a massive galaxy, the source of 
the ionizing radiation, and when or if the gas within the dwarf galaxies is ionized. 

Previous investigation by Kepner et al. (1997) suggest that the external UV radiation in early epochs 
can fully ionize residual gas in typical-mass dwarf galaxies, and so rapid mass loss might be expected. The 
results in §3 show, however, that hot ionized gas can be retained in potentials with velocity dispersions 
comparable to or less than the sound speed of the gas, provided that the host galaxy is isolated. Thus, even 
in the presence of an intense UV background, low-mass dwarf galaxies do not automatically lose their gas 
content, albeit the density distribution can become relatively flat for systems which contain hot gas. 
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Such an environmental dependence may help to explain differences between some types of dwarf galaxies. 
Dwarf irregular (dirr) and dwarf elliptical (dE) galaxies are observed to be similar in terms of their dark- 
matter potentials and ordinary matter content. But, a significant amount of ordinary matter in dirrs is 
in the form of gas, whereas little gas is found in dEs. It has been suggested that this difference is due to 
the efficiency of gas removal (Lin & Faber 1983), motivated by the fact that dEs arc primarily located in 
clusters of galaxies, whereas dIrrs arc mostly found in the field. Although the presence of hot gas in clusters 
of galaxies can lead to efficient stripping (Murray et al. 1993), the present results indicate that tidal effects 
can also remove the ionized gas in dwarf galaxies within a galactic cluster environment within several Gyr. 

At the center of large clusters, however, where the density is nearly homogeneous or declines slowly 

with radius, the tidal potential is compressive rather than disruptive. In these regions, dwarf galaxies may 
also be able to retain their gas content, even if it is mostly ionized. Nucleated dwarfs are often found in 
these central regions of clusters of galaxies. The retention of residual gas may enable the dwarf systems 
to undergo repeated star formation, possibly leading to the formation of their nucleated structure. This 
dependence upon location is in the same sense as that proposed by Babul & Rees (1992). In that earlier 
work, it was proposed that the added confinement by the denser gas near cluster centers would help to 
prevent supernova-driven mass loss from dwarf galaxies, allowing them to retain their gas for subsequent 
generations of star formation. Away from the cluster center, however, where tidal forces are disruptive, 
they shall act in opposition to external confinement, to drive mass loss. In those regions, the presence of a 
significant external medium may even enhance mass loss, via the Kelvin-Helmholtz instability (Murray et 
al. 1993). 

The evolution of non-isolated dwarf systems depends upon both their masses and their proximity to a 
massive galaxy. The results in § 4 indicate that if a massive galaxy forms near a dSph prior to the efficient 
conversion of the internal gas of the dSph into stars, then the remaining gas shall be rapidly lost to tidal 
stripping if it becomes ionized. In the most extreme case, when background UV heating has ionized the gas 
within the dSph, preventing it from forming any stars, tidal stripping of the gas shall lead to the formation of 
a starless, dark-matter mini-halo orbiting the larger galaxy. Only the gas component of the smaller systems is 
efficiently lost, because the outflow is continually replenished by the attempt to re-establish quasi-hydrostatic 
equilibrium. The dark-matter mini-haloes may be preserved, because they were formed in relatively high 
density environments (Moore et al. 1999). Thus, the results presented here naturally imply that the halo 
of large galaxies may be populated with many 'stripped' dark matter, self-gravitating systems. 

By contrast, the gas within the most massive dSphs, and that within dEs may be self-shielded, and 
remain mostly neutral. That, coupled with their deeper potentials, allows them to retain their gas for a 
much longer time than low mass dSphs, and form multiple generations of stars before losing their gas, even 
when in close proximity to a massive galaxy. The dependence of the rate of tidal stripping upon the mass 

of the dwarf system may provide an explanation for the "missing satellite galaxy" problem, because, in 
accordance with the CDM model, the galactic mass function decreases with Md- 

1/2 

Tidal loss is also crucial for globular clusters. As in dSphs, (p^ for globular clusters is a few to a few 
tens of km s^^, comparable to the sound speed of ionized gas. The highly homogeneous chemical abundance 
among stars within any given globular cluster (Sandage & Katem 1977; Cohen 1979; Sandage & Katem 
1982; Richer & Fahlman 1984; Boltc 1987a, b) reqiurcs thorough mixing of the proto-cluster clouds. Even 
with efficient turbulent mixing, many internal crossing time scales are required (Murray & Lin 1990). Thus, 
prior to the formation of presently observable stars in globular clusters, gas must be retained for a time 
comparable to or longer than their Galactic orbital period. The baryonic density of globular clusters is 
higher than that of dSphs, and so, prior to efficient star formation, the gas was likely to be self-shielded from 
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external heating. Because the value of ai is much larger for cool atomic and molecular gas, the gas can then 
be retained for several Galactic orbits. Self-gravity of the gas was also probably very important for globular 
clusters at early epochs, which would further help them to retain their gas. 

The internal UV production following efficient star formation within a globular cluster would, however, 
ionize the gas, after which would be rapidly lost to the system. The same fate would apply to the most 
massive dSphs orbiting close to a massive galaxy. In the presence of a tidal field, therefore, star formation 
within dwarf systems can become a self-terminating process, even in the absence of ram pressure stripping, 
supernova heating, or other means of violent gas ejection. 

In the hierarchical galaxy formation scenario, large galaxies are formed through coalescence of dwarf 
galaxies. Galactic tides enable the rapid decoupling of the gas from dwarf galaxies. The tidally removed gas 
is detached gently from the dwarf systems, and shall form thin tidal streamers, similar to the Magellanic 
Stream. These gas streams are dispersed by differential gravity, precession, nonlinear shock dissipation, and 
drag by the residual gas in the halo. The tidally removed gas does, however, retain its angular momentum. 
When it is mixed with the residual halo gas, together they can form the disk (Sommer-Larsen, Gelato, & 
Vedcl 1999). The "tidally stripped" dark matter clumps shall themselves undergo subsequent orbital decay 
due to dynamical friction. In outer regions of the of the large parent galaxies, the residual gas in the halo is 
tenuous, possibly allowing the streamers to preserve their integrity. Repeated Galactic encounters may lead 
to recapture of the gas from the streamers by the dwarf system, and induce secondary star formation. The 
follow-up investigation on the tidally removed gas will be analyzed and presented in subsequent work. 

This work was performed under the auspices of the U.S. Department of Energy by University of Cal- 
ifornia, Lawrence Livermore National Laboratory under Contract W-7405-Eng-48. This work is partially 
supported by NASA through an astrophysical theory grant NAG5-12151. 
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Table 1. One-Dimensional Galaxy Models 
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Table 2. Three-Dimensional Galaxy Models 



Model 
units 
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Rc 
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1 
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2 
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3 
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Table 3. Three-Dimensional Cluster Models 
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Fig. 2. — Mass of gas contained within the tidal radius as a function of time for the one-dimensional models 
lD-1, lD-2, and lD-3. 
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Fig. 3. — Mass of gas contained within the tidal radius as a function of time for the one-dimensional models 
lD-1, lD-4, and lD-5. 
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Fig. 4. — Density of gas for the one-dimensfonal model. The thicker Unes denote the density of neutral gas 
Pn, while the thinner ones denote the density of ionized gas p,. 
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Fig. 5. — Surface plot of 4>{x,y,z = 0) for the three-dimensional galaxy Model 3D-ld. In code units, the 
y-axis runs from to 1, the x-axis from -1 to 1, and Xt = 0.89. The potential, 0, is in code units. 
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Fig. 6. — The values of M2 (t) /M2 (0) for the three-dimensional dwarf galaxy models. The results for Model 1 
are indicated by the solid curve, Model 2 by the long dashed curve, Model 3 by the short dashed curve, and 
Model 4 by the dot-dashed curve. The values of tioss are found from exponential fits to the curves. 
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Fig. 8. — The values of M2(t)/M2(0) for the cluster models, displayed as in Figure 6. 



